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ABSTRACT 

Recent far-infrared (IR) observations of supernova remnants (SNRs) have revealed significantly large amounts of 
newly condensed dust in their ejecta, comparable to the total mass of available refractory elements. The dust masses 
derived from these observations assume that all the grains of a given species radiate at the same temperature, 
regardless of the dust heating mechanism or grain radius. In this paper, we derive the dust mass in the ejecta of the 
Crab Nebula, using a physical model for the heating and radiation from the dust. We adopt a power-law distribution 
of grain sizes and two different dust compositions (silicates and amorphous carbon), and calculate the heating 
rate of each dust grain by the radiation from the pulsar wind nebula. We find that the grains attain a continuous 
range of temperatures, depending on their size and composition. The total mass derived from the best-fit models 
to the observed IR spectrum is 0.019-0.13 M 0 , depending on the assumed grain composition. We find that the 
power-law size distribution of dust grains is characterized by a power-law index of 3. 5-4.0 and a maximum grain 
size larger than 0. 1 /xm. The grain sizes and composition are consistent with what is expected for dust grains formed 
in a Type IIP supernova (SN). Our derived dust mass is at least a factor of two less than the mass reported in previous 
studies of the Crab Nebula that assumed more simplified two-temperature models. These models also require a 
larger mass of refractory elements to be locked up in dust than was likely available in the ejecta. The results of this 
study show that a physical model resulting in a realistic distribution of dust temperatures can constrain the dust 
properties and affect the derived dust masses. Our study may also have important implications for deriving grain 
properties and mass estimates in other SNRs and for the ultimate question of whether SNe are major sources of 
dust in the Galactic interstellar medium and in external galaxies. 
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1. INTRODUCTION 

In recent years, the far-infrared (IR) observations of super- 
nova remnants (SNRs) with the Herschel Space Observatory 
(Pilbratt et al. 2010) have allowed us to finally tackle the long- 
standing question of whether supemovae (SNe) contribute a 
significant amount of dust to the interstellar medium (ISM). 
Theoretical dust condensation models predict that 0. 1-0.7 M© 
of dust should form in a SN explosion of a ~20M© star 
(Todini & Ferrara 2001; Cherchneff & Dwek 2010; Nozawa 
et al. 2010; Silvia et al. 2012). A significant fraction of this 
newly condensed dust may be destroyed following the encounter 
with the reverse shock in the ejecta (e.g., Dwek 2005; Bianchi 
& Ferrara 2005; Kozasa et al. 2009), and SN blast waves in 
the ISM (Jones & Nuth 2011 and references therein). SNe may 
therefore be required to account for the mass of dust observed 
in local (e.g., Dwek & Scalo 1980; Matsuura et al. 2009; Calura 
et al. 2008; Boyer et al. 2012; Zhukovska & Henning 2013) 
and high-redshift galaxies (e.g., Dunne et al. 2003; Morgan & 
Edmunds 2003; Dwek et al. 2007 ; Michalowski et al. 2011; Gall 
et al. 2011a, 2011b; Dwek & Cherchneff 2011; Valiante et al. 
2012). Required SN dust yields may be as low as ~0. 1 M© or as 
high as ~ 1-2 M©, depending on the grain destruction efficiency 
in the ISM (Dwek et al. 2007). 

Spitzer Space Telescope observations of SNRs have revealed 
ejecta dust in many remnants with estimated masses in the 
0.02-0. 1 M q range (Sugerman et al. 2006; Rho et al. 2008, 2009; 
Temim et al. 2010, 2012b), but the recent far-IR observations 


that are sensitive to cooler dust are now suggesting even larger 
masses. Cas A appears to contain 0.075 M© of cool (~35 K) 
ejecta dust located on interior of the reverse shock (Barlow 
et al. 2010; Sibthorpe et al. 2010), raising the total dust 
mass to 0.1 M© (Rho et al. 2009). Based on recent Herschel 
observations, Matsuura et al. (201 1) reported 0.4-0.7 M© of cool 
dust in SN 1987 A, while Gomez et al. (2012) find 0. 12-0.25 M© 
of dust in the Crab Nebula. 

While these recent results imply that SNe may indeed be 
important suppliers of dust, the masses in each of these cases 
are derived using simple modified blackbody distributions with 
one or two dust temperature components to fit the IR and 
submillimeter spectral energy distribution. A more realistic 
scenario is to identify the heating source in the ejecta, and use 
a continuous size distribution of grains with each grain heated 
to a different temperature that depends on its composition, size 
and optical properties. For example, Richardson et al. (2013) 
modeled the H 2 emission in the Crab Nebula, and showed that a 
given grain composition and size distribution leads to continuous 
distribution of dust temperatures. 

In this work, we present a detailed model of the physical 
mechanism giving rise to the observed IR emission from the 
Crab Nebula. Dust in the Crab Nebula was discovered as an 
IR excess above the synchrotron power-law spectrum of its 
pulsar wind nebula (PWN; Trimble 1977; Glaccum et al. 1982; 
Marsden et al. 1984; Douvion et al. 2001; Green et al. 2004; 
Temim et al. 2006). Absorption features from dust are observed 
to spatially coincide with the ejecta filaments (Woltjer & 
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Veron-Cetty 1987; Fesen & Blair 1990; Hester et al. 1990; Blair 
et al. 1997; Loll 2010), and a recent analysis of the Spitzer 
Infrared Spectrograph (IRS) spectra confirmed that the dust is 
indeed located in the filament cores (Temim et al. 2012b), and 
has therefore condensed in the ejecta. 

The model compares three different dust compositions, each 
characterized by a power-law distribution in grain radii whose 
index is allowed to vary. The primary goal of our study is 
to determine whether a physical model for the IR emission, 
characterized by a continuous distribution in dust temperatures, 
will affect the derived grain properties and dust mass in this 
remnant. Our choice of the Crab Nebula is primarily motivated 
by the fact that the dust grains are heated by the synchrotron 
radiation from the PWN, so that their temperature can be 
accurately derived for each size and composition. Furthermore, 
since the newly formed dust has neither been processed by the 
reverse shock, nor mixed with the ambient medium (for a review, 
see Hester 2008; Loll 2010), the Crab Nebula offers a unique 
opportunity to study the mass, composition, and size distribution 
of pristine SN-condensed dust. 

In spite of the fact that we concentrate on modeling the IR 
emission from the Crab Nebula, our results may have broader 
implications for deriving dust properties and masses from IR 
observations of SNRs. 


2. DUST HEATING MODEL 
2.7. Heating by the PWN 

Since the dust in the Crab Nebula is concentrated in the 
Rayleigh-Taylor filaments that form a cage around the PWN 
(Temim et al. 2012b), we model the heating of dust by assuming 
that the heating source is located at the center of the Crab Nebula, 
at the location of the Crab pulsar. The dominant heating source 
for the dust in the Crab Nebula is the synchrotron radiation 
from the PWN, with an insignificant possible contribution from 
collisional heating by the gas in the filaments (Dwek & Werner 
1981 ; Temim et al. 2012b). In this case, the heating rate (erg s -1 ) 
of a single dust grain of radius a is given by 


H(a ) = 


7ra 2 / U(v)Q abs (v, a)dv 
4nr 2 


( 1 ) 


where L v is the non-thermal specific luminosity of the Crab 
Nebula’s PWN, a is the grain size, ga bs(u a) is the absorption 
coefficient for a given grain composition, and r is the distance 
between the radiation source and the dust grain. The non-thermal 
luminosity L v (y) of the Crab Nebula’s PWN that was used in our 
model is summarized in Figure 2 of Hester (2008). As described 
in Temim et al. (2012b), we consider L v {y) up to an energy of 
about 0.6 keV, since we find that the fraction of the energy 
deposited in the dust at higher energies makes a negligible 
contribution of less than 1%, due to the combined effects of 
decreasing of both L v {y) with energy, and the efficiency of the 
energy deposited by the photoelectrons in the dust grain (Dwek 
& Smith 1996). 

Based on the three-dimensional models of the Crab Nebula, 
the ejecta filaments are located between 0.55-1.0 pc from the 
center of the nebula (e.g., Cadez et al. 2004). In our dust model, 
we allow the distance from the heating source r to vary between 
0.5-1 .5 pc in intervals of 0.2 pc. We allow the upper limit on the 
distance to extend beyond the physical location of the filaments 
in order to account for any attenuation by dust that would affect 
the heating rate. The best-fit models described in Section 3 favor 
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Figure 1 . Mass absorption coefficients as a function of wavelength where 
the silicate grains (Weingartner & Draine 2001) are shown in red, amorphous 
carbon- AC grains (Rouleau & Martin 1991) in blue, and amorphous carbon-BE 
grains (Zubko et al. 2004) in green. 

(A color version of this figure is available in the online journal.) 


the low end of the distance range, suggesting that the internal 
absorption in the nebula is negligible. In Temim et al. (2012b), 
we find that the optical depth is indeed low, with r < 1 . 

In equilibrium, the heating rate in Equation (1) of each grain 
of radius a is equal to the radiative cooling rate given by 

L gr (a) = 4 7ta 2 J nB v (v, T)Q ahs (v, a)dv , (2) 

= 4 md(a) J ttB v (v, T)i<(v,a)dv, 

where B v is the Planck function, T is the dust grain temperature, 
k(v , a) = 3 Qabs/4pa is the mass absorption coefficient, = 
4npa 3 /3 is the dust mass, and p is the mass density of a dust 
grain. At the wavelengths of the IR emission where A a (the 
Rayleigh regime), k is independent of grain radius. We used this 
relationship to compute the temperature for each grain radius 
and composition, and each distance r, and produced a final set 
of dust emission models by convolving the resulting emission 
spectrum for each grain size with the size distributions functions 
described in the previous section. 

2.2. Dust Composition 

The mid-IR spectrum of dust in the Crab Nebula is fairly 
featureless (Temim et al. 2012b), consistent with the generally 
featureless spectra generated by silicate and carbon grains. 
Theoretical dust condensation models do indeed show that a 
large fraction of dust formed in Type Up SNe is in the form of 
silicate and carbon (e.g., Kozasa et al. 2009). In our models, 
we use three different sets of optical constants to calculate the 
mass absorption coefficients: those from Weingartner & Draine 
(2001) to characterize the silicate grains; those from Rouleau 
& Martin (1991) to characterize the amorphous carbon grains 
(labeled AC throughout the paper); and those from Zubko et al. 
(2004) for amorphous carbon grains of type BE. The values 
of a: as a function of wavelength for each of the three grain 
compositions is shown in Figure 1 . 

The real and imaginary parts of the refractory index of AC 
carbon, n and k , extend only to 300 /rm. They smoothly increase 
in the 10-300 /rm wavelength range, following the functional 
form that can be approximated by a power law A -0 06 and A -015 
for n and k , respectively. These power law indices were used 
to extrapolate the optical constants out to 10 4 /rm, leading to a 
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Figure 2. Best-fit dust grain heating models for three different grain compo- 
sitions (see Figure 1). The data include the average Spitzer IRS dust spectrum 
(Temim et al. 2012b), and the dust fluxes from Herschel PACS and SPIRE 
(Gomez et al. 2012). The line-free regions of the IRS spectrum that were used 
in the fit are shown as black data points. The best-fit parameters are summarized 
in Table 1, and corresponding total dust masses in Table 2. 

(A color version of this figure is available in the online journal.) 


smooth k which maintains a constant emissivity index at long 
wavelengths (see Figure 1). 

While the optical constants for carbon-BE area measured 
beyond 500 /tm, there are no data below ~0.1/xm, where 
a significant fraction of the energy absorption takes place. 
Therefore, dust temperatures computed for carbon-BE grains of 
all sizes will be lower than the physical temperatures expected 
from PWN heating. We included them in this work only for 
comparison with previous work by Temim et al. (2012b) and 
Gomez et al. (2012). Because of the lack of short wavelength 
coverage for the optical constants of carbon-BE, we chose 
carbon- AC to represent the carbon dust in the ejecta. 

2.3. Grain Size Distributions 

The determination of the grain size distribution in SNRs 
requires detailed knowledge of the physical condition in the SN 
ejecta, and the dust heating mechanism. The size distribution 
of SN condensed dust is therefore a priori unknown. Models 
for the formation of dust in SNe (Todini & Ferrara 2001; 
Nozawa et al. 2010) have derived the size distribution of the 
various dust species that form in SN ejecta. However, these 
calculations assumed a uniform ejecta, and that the growth of 
the dust grains proceeded only through the accretion of single 
monomers, leaving out possible growth through coagulation in 
the ejecta. 

Considering the above-mentioned uncertainties, we adopt 
a general parameterization for the grain size distribution, de- 
scribed by a power law in grain radii, a~ a , and a lower and up- 
per grain radius cutoff on the grain radii, a m i n and a mdLX , respec- 
tively. Historically, such characterization was used by Mathis 
et al. (1977, MRN) to derive the size distribution of interstellar 
dust from the observed UV-optical extinction (for a review see 
Clayton et al. 2003). The incorporation of additional observa- 
tional constraints: the diffuse IR emission, interstellar abun- 
dances, and diffuse interstellar scattering, have yielded a more 
complex interstellar grain size distribution (Weingartner & 
Draine 2001; Li & Draine 2001; Zubko et al. 2004). This dis- 
tribution is determined by the size distribution of the grains that 
are injected from all sources into the ISM, and by the various 
interstellar processes that alter their sizes, including thermo- 
kinetic sputtering, evaporative and shattering grain-grain 


Table 1 

Best-fit Model Parameters 


Composition 

r 

(pc) 

a 

^max 

(M m) 

x 2 

Silicates 

0.5 (0.5) 

3.5 (3.5) 

5.0 (>0.6) 

3.13 

Carbon-AC 

0.7 (0.5-0.7) 

4.0 (3. 5-4.0) 

0.1 (>0.1) 

1.86 

Carbon-BE 

0.5 (0.5) 

3.5 (3.5) 

0.6 (>0.3) 

3.14 


Notes. Best-fit model parameters where r is the distance of the dust grains from 
the center of the PWN, a is the power-law index on the grain size distribution, and 
a max is the maximum grain size cut-off. The values in parentheses represent the 
range of parameter values obtained from the bootstrap method (see Section 2.5). 


collisions, accretion and coagulation. Attempts to character- 
ize the net size distribution resulting from all these processes 
were made by Liffman & Clayton (1989), O’Donnell & Mathis 
(1997), and Hirashita (2012). 

In order to fit the mid and far-IR dust emission in the Crab 
Nebula, we constructed a grid of grain size distributions with the 
power-law index a ranging from 0. 0-6.0, and < 2 max ranging from 
0.03 to 5.0 /xm. The maximum limit of 5.0 /xm is already larger 
than the maximum grain size obtained from dust condensation 
models (e.g., Todini & Ferrara 2001), and grain sizes larger than 
this are not expected to form in SN ejecta. By allowing a max to 
vary, we hope to determine to what radius SN grains can grow, 
and to compare the best-fit parameters to theoretical models for 
grain growth in Type IIP SNe. 

Fitting the minimum grain size cut-off a m { n allows us to 
estimate the amount of future grain processes in SNRs. However, 
in the Crab Nebula the reverse shock has not yet reached the 
PWN, so we do not expect that any grain destruction has yet 
occurred. We tested the effect of varying the minimum grain 
size on the best-fit parameters in our model by allowing a^n to 
be as high as 0.03 /xm the lowest allowed value of a max . Due 
to the lack of short-IR dust spectra, the a m i n parameter is not 
well constrained by our model. However, we did confirm that 
varying this parameter does not affect the best fit dust mass and 
shape of the grain size distribution, since most of the dust mass 
is contained in the larger grains. We therefore fixed the value of 
<2 m in 1° be 0.001 /xm. 

2.4. Fitting of the IR Spectrum 

We fitted our entire grid of models to the observed dust emis- 
sion from the Crab Nebula that include the average Spitzer IRS 
spectrum of the dust emission scaled to the total synchrotron 
and line- subtracted MIPS 24 /xm flux (Temim et al. 2012b), and 
the synchrotron and line- subtracted Herschel PACS 70, 100, 
and 160 /xm flux measurements, SPIRE 250 and 350 /xm mea- 
surements, and revised Spitzer fluxes from Gomez et al. (2012), 
where cold dust emission still contributes. Even though the 
line emission was subtracted from the entire IRS spectrum, our 
fits only included the wavelength regions of the IRS spectrum 
where the line emission did not contribute significantly. The data 
and the best-fit model for each grain composition are shown in 
Figure 2. The portions of the IRS spectrum data were used in 
the fit are overplotted in black. 

The best-fit values for the distance from the heating source 
r, power-law size distribution index a , and the maximum grain 
size cut-off < 2 max are summarized in Table 1. The dust model 
for the amorphous carbon- AC grain composition provides the 
best fit to the data, since it simultaneously provides a good fit 
to the mid-IR spectrum and the Herschel data points at far-IR 
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Figure 3. The range of spectra simulated by the Monte Carlo bootstrap method 
described in Section 2.5 is shown as the gray band, with the actual data overlaid 
in black. 


wavelengths. The best-fit silicate grain composition falls some- 
what short at far-IR wavelengths, indicating that the presence 
of colder grains would be required to fit the far-IR data points. 
This would either require a grain radius larger than our limit 
of 5 /x m, or it may mean that some of the dust may be located 
farther away from the heating source than the best-fit average 
distance. 

2.5. Monte Carlo Simulations 

We note that the large uncertainties in the data points in 
the IRS spectrum are dominated by systematic uncertainties 
introduced by the subtraction of the underlying synchrotron 
spectrum (see Temim et al. 2012b). For this reason, the reduced 
X 2 for our best fits is too small, on the order of ~0.06. Since our 
uncertainties are not random statistical uncertainties, we cannot 
use x 2 statistics in determining confidence levels for our fitted 
parameters. Instead, we used the bootstrap method to estimate 
the spread in the fitted parameters. The absolute x 2 values from 
our fits are used only to demonstrate the relative goodness of fit 
for the various models. 

In order to estimate the uncertainties in the fitted parameters 
while correctly accounting for the statistical uncertainties in the 
dust spectrum introduced by the subtraction of the synchrotron 
spectrum, we generated 5000 simulated spectra using the 
bootstrap method. The spectra were simulated by adding an 
average global synchrotron spectrum to select IRS data (black 
data points in Figure 2), with a synchrotron power-law index 
of 0.42, as found by Gomez et al. (2012). We then produced 
5000 sample dust spectra by resubtracting a global synchrotron 
spectrum while randomly selecting spectral indices within 3 a 
(±0.02) of the best fit value of 0.42. We also randomly selected 
values for the 3.6 jam flux from which the synchrotron spectrum 
is extrapolated in the range of 12.65 ±0.63 Jy (Temim et al. 
2012b), and we accounted for the ~5% uncertainty that results 
from normalizing the IRS dust spectrum to the total flux from 
dust in the MIPS 24 jam band. The Spitzer MIPS and Herschel 
data were similarly varied within the uncertainties shown in 
Figure 4 of Gomez et al. (2012). The resulting range of simulated 
spectra are shown as the gray band in Figure 3, with the actual 
data overlaid in black. We then fitted our grid of dust models to 
each of the simulated spectra to obtain a distribution in the best- 
fit values for the power-law index a of the grain size distribution, 
the value of a max , and the total dust mass for each composition. 
The results for each parameter will be discussed in the following 
sections. 



0.001 0.010 0.100 1.000 10.000 
Grain Size (/xm) 

Figure 4. Grain temperature as a function of grain size for radiatively heated 
dust in the Crab Nebula, located 0.5 pc from the center of the PWN. The input 
spectrum of the heating source is the broadband non-thermal spectrum of the 
Crab Nebula, as summarized in Figure 2 of Hester (2008). The colored curves 
correspond to grains with mass absorption coefficients shown in Figure 1 . As 
described in Section 3.1, the large difference in the temperatures of the two 
carbon compositions is primarily due to the lack of short-wavelength data for 
the absorptions coefficients of BE carbon. 

(A color version of this figure is available in the online journal.) 

3. RESULTS AND DISCUSSION 
3.1. Distribution of Grain Temperatures 

The main difference between using a physical dust heating 
model and a simple two-temperature fit to the data is that in 
the former, the range of dust temperatures and their distribution 
is uniquely determined by the grain size distribution and the 
heating mechanism. Figure 4 shows the dust temperature (7) as 
a function of grain size, assuming that the heating source (i.e., 
the Crab’s PWN) is located at the best-fit distance of 0.5 pc. 
The colored curves show the equilibrium grain temperature as 
a function of grain size for each of the three grain compositions 
that were used in the fit. The shape of the temperature profile 
was obtained by equating the dust heating and cooling rates, as 
described in Section 2.1. It is important to note the significant 
difference in the temperature of carbon-AC and carbon-BE 
grains. The additional coverage between 0.002 and 0.05 /xm 
(0.02-0.6 keV) for the carbon-AC grains raises the temperature 
for all grain sizes, but has a most apparent affect on the 
smaller grains (see Figure 4). While the carbon-AC grains with 
sizes below ~0.01 /xm are heated to a temperature of ~68 K, 
the same- sized carbon-BE grains only attain a temperature of 
~55 K. The large difference in dust temperatures between these 
two carbon compositions is not due to the differences in grain 
properties or the shape of the absorption coefficients gabs> but 
is instead primarily due to the limited wavelength coverage of 
gabs- A significant fraction of the radiative energy of the PWN 
is emitted at energies above 0.02 keV (0.05 /xm), where there 
are no data for the absorption efficiencies of the carbon-BE dust. 
Accounting for energies above 0.02 keV is especially important 
for heating by PWNe, since their non-thermal spectra peak at 
these energies. So even if we cannot determine the nature of the 
carbon dust in the Crab, whether it is AC or BE, our choice of 
carbon-AC dust is primarily motivated by the fact that it gives a 
more physical representation of the temperature distribution in 
the nebula. 

Figure 5 shows the comparison between a two-temperature 
fit to the dust emission in the Crab Nebula using a silicate 
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Figure 5. Top: best-fit two-temperature fit for the silicate grain composition. 
The blue curve represents the warm grains with a temperature of 55 ± 11 K, 
while the red curve represents the cold 30 ± 5 K grains. The black curve is the 
sum of the individual components. Middle: our best-fit dust heating model for 
silicate grains. The rainbow colored curves represent the spectra of individual 
dust grains of different sizes that are heated to different temperatures. The black 
curve represents the sum of the individual rainbow-colored spectra. Bottom: 
the same model as in the middle plot, but with the individual spectra added 
cumulatively starting from the hottest grains. 

(A color version of this figure is available in the online journal.) 


Figure 6. Top: best-fit one-temperature fit for the carbon- AC grain composition. 
The blue curve represents a grain temperature of 56 ± 2 K. Middle: our best-fit 
dust heating model for carbon- AC grains. The rainbow colored curves represent 
the spectra of individual dust grains of different sizes that are heated to different 
temperatures. The black curve represents the sum of the individual rainbow- 
colored spectra. Bottom: the same model as in the middle plot, but with the 
individual spectra added cumulatively starting from the hottest grains. 

(A color version of this figure is available in the online journal.) 


grain composition (top panel), and our best-fit physical dust 
heating model with parameters listed in Table 1 (middle and 
bottom panels). The middle panel of Figure 5 shows spectra of 
grains with different sizes as individual rainbow colored curves 
that all add up to the best-fit black curve. The bottom panel 
shows the same model, but with every tenth cumulative sum of 


individual spectra shown in rainbow colors. While both models 
provide equally good fits to the IR data, the physical dust heating 
model allows us to constrain physical properties of the dust, and 
reduces the total dust mass responsible for the IR emission (see 
Section 3.3). Figure 6 shows the same three plots, but for the 
carbon- AC grain composition. We note that the IR spectrum 
can be well fitted with only a single temperature component 
of carbon- AC. The physical model still produces a range of 
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Figure 7. x 2 contour plots for the maximum grain size cutoff, a max , and the 
index on the power-law grain size distribution a. The contours indicate that an 
a value of 3. 5-4.0 and a larger grain size cutoff tend to produce better fits for 
both silicate and carbon- AC grain compositions. 

(A color version of this figure is available in the online journal.) 

temperatures, but from the comparison of Figures 5 and 6, it 
is evident that the best fit physical model produces a narrower 
temperature range for carbon- AC than for silicates. 

3.2. Grain Size Distribution 

The dust heating model also allows us to place some con- 
strains on the shape of the grain size distribution and the max- 
imum size of grains that formed in the Crab Nebula. Figure 7 
shows the x 2 contour plots for a , the index of the power-law 
size distribution, versus the maximum grain size cut-off a m ax . 
The parameter values for models with minimum x 2 values are 
listed in Table 1. The contour plots show the relative good- 
ness of fit for the range of parameter values used in our grid 
of dust distribution models. For all compositions, the best-fit 
value for the power-law index is a = 3. 5-4.0. The fits to the 
5000 Monte Carlo simulations show that the power-law index 
is very well constrained with a standard deviation of only 0.1. 
This size distribution is similar to what is found for the ISM 
(e.g., Mathis et al. 1977). Physical modeling of dust emission 
for a larger sample of SNRs is needed to determine if similar 
size distributions are found in other core-collapse SNe. 

In Temim et al. (2012b), we found that the dust emission 
in the Crab observed by Spitzer implies a small grain size of 
<0.05 /z m. However, since we only had coverage up to 70 /z m, 
we were sampling only the smaller grains that were heated to 
temperatures of around 50 K (see Figure 4). The addition of 
the Herschel long wavelength data, and a better constrains on 
the synchrotron index derived from the Planck data (Gomez 


Composition 


Models 

Nucleosynthesis Yields 


Physical 

One/Two-temp 

WW95 

N06 

WH07 

Silicates 

0.13 ±0.01 

0.2 ± 0.1 

0.08 

0.32 

0.12 

Carbon-AC 

0 019 +001 ° 
U.U 1 V — 0003 

0.018 ±0.005 

0.05 

0.10 

0.11 

Carbon-BE 

0 040 + ° 021 
u.m-u_ 00 io 

0.08 ± 0.03 





Notes. The uncertainties on all dust masses in the table were determined from 
bootstrap Monte Carlo simulations of the data (see Section 3.3). The dust masses 
listed for the nucleosynthesis models are the maximum allowed values assuming 
a 100% condensation efficiency in the ejecta, and the yields for a 11 Mq 
progenitor from Woosley & Weaver (1995, WW95), and a 13 Mq progenitor 
from Nomoto et al. (2006, N06) and Woosley & Heger (2007, WH07). 

et al. 2012), has now allowed us to place some constraints on 
the maximum size of dust grains in the Crab’s filaments. Our 
best-fit model appears to favor a fairly large maximum grain 
size cut-off (see Figure 7). The value for a max for the best-fit 
model is 0.1 /z m or larger, depending on the grain composition 
(Table 1). This is consistent with models of dust formation 
in Type IIP SNe (Kozasa et al. 2009; Nozawa et al. 2010). 
These models suggest that the dust mass of grains formed in 
Type IIP SNe with massive H-envelopes is dominated by grains 
larger than 0.03 /zm, while the mass of dust for Type lib SNe is 
dominated by very small grain <0.006 /zm. Probing the grain 
size in SNRs through physical heating models can therefore be 
a useful tool for comparing and constraining dust condensation 
models. The Crab Nebula is thought to be a result of a Type IIP 
explosion, primarily based on the low expansion velocities of 
the ejecta filaments and the low progenitor mass (e.g., Chevalier 
2005; MacAlpine & Satterfield 2008). The large dust grain radii 
inferred from our models are consistent with a Type IIP origin. 
We note that the larger grain size also implies that a large 
fraction of the dust mass will survive the eventual reverse shock 
interaction (Dwek 2005; Kozasa et al. 2009; Bianchi et al. 2009; 
Nozawa et al. 2010). 

3.3. Revised Dust Mass for the Crab Nebula 

In order to compare our derived dust masses to previous 
studies, we first reproduced the two-temperature fits from 
Gomez et al. (2012) with the same dust grain compositions 
that they used in their analysis. Our best-fit temperatures are 
consistent with the values of Gomez et al. (2012). The total 
masses that we derive are 0.2 zb 0.1 M 0 and 0.08 d= 0.03 M 0 
for silicate and carbon-BE grains, respectively. They are listed 
in Table 2, and also found to be consistent with the values 
from Gomez et al. (2012) of 0.25 +q 3 0 2 8 M q for silicates and 
0.12 ± 0.01 Mq for carbon-BE grains. 

We then calculated the total dust mass corresponding to our 
best-fit dust heating models characterized by a grain size and 
temperature distribution, the latter obtained by calculating the 
temperature of each dust grain as it is radiatively heated by the 
PWN. We find that the total mass of the silicate and carbon- 
BE dust is reduced by a factor of ~2. We derive a dust mass 
of 0.13 ± 0.01 Mq for silicates and 0.04+ f° r carbon- 
BE grains. The total mass for the carbon-BE model is presented 
only to demonstrate the effect of using a continuous temperature 
distribution, as opposed to a bimodal temperature distribution, 
on the calculated dust mass. As discussed in previous sections, 
the lack of short-wavelength data for optical constants for 
carbon-BE grains prevents us from estimating the actual grain 
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Figure 8. Distribution of total dust masses obtained from the bootstrap fits to the 
5000 simulated spectra using the silicate (top) and carbon- AC grains (bottom). 
The dotted vertical lines represent limits for which 90% of masses lie either 
above or below the corresponding mass. These mass limits are taken as the 
uncertainty range in our best-fit total masses (see Table 2). 



temperature for this composition. Consequently, the resulting 
dust properties are not a physical representation of their actual 
values in the Crab Nebula. For this reason, we chose to model the 
carbon emission using carbon- AC dust, which produces a more 
realistic distribution of dust temperatures for the PWN-heated 
dust. 

The total dust mass derived from our model for carbon- AC 
is 0.019 + qixb M o- The difference in the mass between the two 
different carbon compositions is in part due to the difference 
in their optical constants at long wavelengths. However, the 
comparison between the masses derived for these two different 
compositions is unphysical, since the lack of optical data for 
the carbon-BE dust at X-ray and UV wavelengths significantly 
affects the derived temperature distribution, and hence the mass. 
Surprisingly, we find that the IR data are well fitted even with a 
single temperature component of carbon- AC dust with the same 
mass as derived from our physical model and a temperature of 
T = 56 it 2 K (see Table 1 and Figure 6). This shows that our 
parameterization of the grain size distribution is flexible enough 
to produce a temperature distribution leading to the same dust 
mass as that derived from a simple single-temperature fit. It 
also suggests that a distribution of grain temperatures does not 
automatically lead to a significant change in dust mass and that, 
for a given IR spectrum, the magnitude of the effect depends on 
the dust composition. 

Table 2 also lists the uncertainties in the total masses that were 
derived using the bootstrap method described in Section 2.5. 
Figure 8 shows the histograms of the total dust mass corre- 
sponding to the best-fit model for each of the 5000 simulated 


spectra. The dotted vertical lines represent limits for which 90% 
of masses lie either above or below the corresponding mass. 
We take the uncertainty range in our derived masses to be the 
range between these two limits, and we list them in Table 2. 
The mass uncertainties in our two-component fits were also 
determined by the bootstrap method. The comparison between 
the mass uncertainties shows that the total mass determined by 
the physical model is better constrained than the simplified one 
and two temperature models, for both silicate and carbon grain 
compositions. Overall, the strict limits placed on the tempera- 
ture distribution by our PWN heating model allow for a more 
physical and accurate determination of the total mass. The ab- 
solute uncertainties on the total mass are likely larger due to the 
uncertainty on the distance to the Crab Nebula, however, this is a 
systematic uncertainty that will equally affect both the physical 
and two-temperature mass estimates. 

Carbon- AC dust provides a somewhat better fit to the IR 
emission in the Crab Nebula than silicate dust, suggesting that 
the dust mass may either be dominated by carbon dust, or may 
possibly be a mixture of the two dust compositions. The best- 
fit silicate model falls short at longer wavelengths and would 
require an unrealistically large grain size (larger than our limit 
of 5 /nm) to reach cold enough temperatures to better fit the 
long wavelength data. A second dust component is therefore 
required to provide a better fit to the data, suggesting that the 
composition is more likely to be a mixture of carbon and silicate 
dust. This is also supported by the fact that the derived mass of 
silicate dust is in excess of the total mass of available refractory 
elements needed to condense these grains. Table 2 lists the 
nucleosynthetic yields from core-collapse SNe for progenitors 
in the 11-13 M 0 range (Woosley & Weaver 1995; Woosley & 
Heger 2007; Nomoto et al. 2006). We assumed a 100% grain 
condensation efficiency in the ejecta, and calculated the upper 
limit on the carbon and silicate dust yield from the total yield 
of carbon, oxygen, magnesium, silicon, and iron. To obtain the 
upper limit on the silicate grain mass, we summed the maximum 
possible masses for Si02, MgO, and Fe 3 C >4 grains. The results 
show that the dust masses derived from the two-temperature 
component models are clearly too large, especially considering 
that the progenitor of the Crab Nebula is estimated to be around 
~9.5 Mq (MacAlpine & Satterfield 2008). The total mass of 
silicate dust derived from the physical models is also too large 
compared to most nucleosynthesis yields, suggesting that the 
dust in the Crab Nebula is predominantly in the form of carbon 
dust. 


4. CONCLUSIONS 

We modeled the mid and far-IR dust emission from the Crab 
Nebula in order to determine the grain properties and total dust 
mass in the ejecta. We use a physical dust heating model in 
which a continuous power-law size distribution of dust grains 
is radiatively heated by the PWN, giving rise to the thermal 
IR continuum. The best-fit models to the dust emission favor 
a grain size distribution with a power-law index of 3. 5-4.0, 
and a relatively large maximum grain size cut-off of >0.1 /rm, 
consistent with theoretical predictions of dust formation in 
Type IIP SNe (Kozasa et al. 2009; Nozawa et al. 2010). 

The IR emission is best described by an amorphous carbon 
grain composition with a total mass of 0.019+ 0003 ^o> or a 
silicate grain composition with a mass of 0.13 zb 0.01 M 0 . The 
dust is likely to be a mixture of carbon and silicate grains, 
with a total mass between these two values. Both values are 
lower than previous estimates from more simplified one- or 
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two-temperature models. The difference in the dust mass is due 
to the use of a continuous grain size and temperature distribution 
derived from a physical model, and a different choice of dust 
compositions with a more complete coverage of their optical 
constants over the range of energies at which the heating by the 
PWN is important. 

The limits placed on the temperature distribution by our phys- 
ical heating model allow for a more accurate determination of 
the total mass. Even though a distribution of grain temperatures 
may not necessarily lead to a significant change in dust mass, as 
is the case for carbon- AC, modeling the data with a physical dis- 
tribution of temperatures is required to determine the magnitude 
of the effect. 

In recent years, most large dust mass estimates for SNRs as- 
sumed the simplified one or two-temperature models, including 
the mass estimates from Herschel data for Cas A, SN 19 87 A, 
and the Crab Nebula (e.g., Barlow et al. 2010; Sibthorpe et al. 
2010; Matsuura et al. 2011; Gomez et al. 2012). Our modeling 
results imply that a physical dust emission model is important 
for deriving a the actual mass of dust in SN ejecta. 

Such models will also provide information on the size 
distribution of the grains, which is important for estimating the 
ultimate survival of the dust following the encounter with the 
reverse shock. The amount of SN-condensed dust in the Crab 
Nebula is significantly lower than that required to produce the 
massive amount of dust observed in high-redshift galaxies even 
if the dust is not destroyed in the ISM. 

The importance of using physical models to determine dust 
properties in SN ejecta is not limited to cases of radiatively 
heated dust. Modeling of collisionally heated dust in, for 
example, SN1987A (Dwek et al. 2010), Puppis A (Arendt et al. 
2010), Kes 75 (Temim et al. 2012a), or LMC remnants (Williams 
et al. 2011), also provide, in addition to dust masses and grain 
size distributions, information on grain sputtering efficiencies 
in the X-ray emitting plasmas, as well as a diagnostics of the 
plasma conditions in the shocked gas. Physical models are thus 
crucial for deriving dust masses and dust properties for other 
PWNe and SNRs. Such studies are important for determining 
the role of SNe as dust sources in the local and high-redshift 
universe. 

We thank Rick Arendt and George Sonneborn for useful 
discussions and comments on the paper, and an anonymous 
referee for comments which led to valuable improvements in 
the paper. 
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